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SUMMARY ^ 

By numerical solution of Busemann's Second Order Airfoil Theory and the 
spanwise summing of airfoil strips, FIN determines the pressure coefficient 
distributed over a given fin configuration moving at supersonic speeds. In de- 
termining the distribution, the program can include the effect of a fin-tip Mach 
cone. 

From this basic calculation, FIN can determine as functions of angle of at- 
tack the lift coefficient, wave drag coefficient, pitching moment coefficient, and 
center-of-pressure location, and as a function of fin cant angle the rolling mo- 
ment coefficient. FIN can also determine the lift coefficient slope, wave drag at 
zero angle of attack, pitching moment coefficient slope, rolling moment coeffi- 
cient slope, and center-of-pressure location at zero angle of attack as functions 
of Mach number. 

Comparisons with windtunnel data show that predicted values using FIN out- 
put fall well within 10% of experiment results. 
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A COMPUTER PROGRAM FOR CALCULATING 
THE AERODYNAMIC CHARACTERISTICS OF 
FINS AT SUPERSONIC SPEEDS 


BACKGROUND 

Calculating the aerodynamic characteristics of the fins of a sounding rocket 
is one of the basic steps in calculating the overall aerodynamics of a sounding 
rocket. Because sounding rockets fly at supersonic speeds for the greater por- 
tion of their flight, the regime at these speeds is of particular interest. 

Methods previously used for calculating supersonic fin aerodynamics 
(references 1,2, 3, and 4) are lacking in accuracy and applicability because of 
the many assumptions and approximations inherent in using them to reach a 
closed form or nearly closed form solution. 

By using a high-speed computer to numerically solve basic theoretical 
equations, one may obtain answers rapidly and as accurately as desired. The 
only restrictions to accuracy are the assumptions inherent in deriving the basic 
equations. The theory found to be most amenable to programming is Busemann’s 
Second Order Supersonic Airfoil Theory as described in reference 1. 

PROGRAM CAPABILITIES AND LIMITATIONS 

Given a supersonic fin with chord plane symmetry, at a given Mach number, 
FIN computes as functions of angle of attack (a) the following: 

• Lift coefficient, c L 

• Wave drag coefficient, r 

• Pitching moment coefficient, c 

m 

• Center-of-pressure coordinate measured from the reference axis, x 

• Center-of-pressure coordinate measured from the root chord, y 

Since a as defined in the program is equivalent to the fin cant angle (8) at a 
= 0, the rolling moment coefficient c { is considered a function of $. 
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FIN also computes as functions of Mach number the following: 

• Lift coefficient slope, C L 

• Wave drag at zero angle of attack, C DW 

o 

• Pitching moment coefficient slope, C_ 

“a 

• Rolling moment coefficient slope, C t 

$ 

• C enter- of -pre s s ur e coordinate measured from the reference axis at 
zero angle of attack, X 0 

• Cente r-of-p res sure coordinate measured from the root chord at zero 
angle of attack, Y 0 

The range of a is from zero to a max , as determined by the user, in increments 
of 1 degree. The program allows for a maximum of 100 Mach number points, 
the user specifying the range and increment of the Mach number. 

Since FIN has been designed for use with a four-finned vehicle, the output 
values of C L , C DW , C m , C L , C DW , and C m are for a pair of identical fins; the 
output values of andC“ are for two pairs of perpendicular fins. The values 
of X , Y , X 0 , and Y o apply only to a single typical fin. 

Busemann's Second Order Airfoil Theory is applied subject to the restrictions 
outlined in reference 1, pages 192 to 241, with special attention to Figures 10.3 
and 10.28. However, the use of the third-order terms in this program slightly 
enlarges the applicability shown in this reference. Busemann's theory is applied 
to small streamwise airfoil strips, and the strips are summed in a spanwise 
direction. In addition, the fin-tip Mach cone is accounted for by applying a cor- 
rection factor to the necessary portion of each strip. This correction-factor 
technique was obtained from references 1 and 3. The fin-tip Mach cone correc- 
tion may or may not be included at the user's discretion. 

Figure 1 shows the types of fins to which FIN is applicable. These configur- 
ations cover all of the shapes normally used on sounding rockets and missiles. 

The program input pattern for each type is given in the section on usage. A 
listing of FIN in FORTRAN IV is given in Appendix A. 

AERODYNAMIC THEORY 

Busemann's Second Order Airfoil Theory has been applied to two-dimensional 
airfoil strips with the following assumptions: 
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• All parts of the airfoil surface are in supersonic flow and make small 
angles with the flow. This implies low angles of attack. 

• The leading edge is sharp. The trailing edge must be sharp for a 
double wedge fin or a modified double wedge. For a single wedge and 
a modified single wedge, the effect of the blunt base is neglected. 

• The shock waves are all attached. 

• Each region of flow over the surface acts independently of the others. 


The basic local pressure coefficient equation of the Busemann theory is a 
third-order series expansion: 




K l 7 ll + K 2 rj? + Kj 77 ? - K* 


( 1 ) 


where 


2 

1 


( 2 a) 


_ (y + 1 )M 4 - 4 /3 2 
4/3 4 


( 2 b) 


(y + 1 )M 8 + (2y 2 - 77 -5 )m 6 + 10(y + 1)M 4 - 12M 2 + 8 (2 c) 

6/S 7 


If the flow over the surface region inclined at 17 j to the flow has been pre- 
ceded by a single compressive shock wave anywhere upstream in the flow, the 
value of K* is: 


K * = (y + 1)M 4 [(5 - 3y)M 4 + 4(y - 3)M 2 + 8] <3 

48 

If the flow upstream of the region in question has contained only expansive 
shock-free flow, K* = 0. rj L is the inclination of the leading edge surface. 
Essentially, the first three terms of equation 1 represent the first-, second-, 
and third-order pressures on the surface region respectively. The K term is 
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the third-order irreversible pressure rise through the bow shock. The index, 
i , indicates the surface for which the pressure coefficient is being calculated. 

The pressure coefficient of equation 1 is applied to a typical airfoil region, 
including the fin-tip Mach cone, by dividing the region into two portions at the 
intersection of the airfoil strip and the Mach cone (Figure 2). A correction factor of 
l A is applied to the local pressure coefficient for the region within the Mach cone. 
The center of pressure of each portion is taken as the midpoint of the length of 
that portion. By this means, equations are obtained for the local lift(F 1 ) , 
drag (F d ) , local center of pressure ( X i ) , and, hence, the local pitching moment 

(«p,) : 1 



The airfoil-strip lift (Fj ) , drag (F d ), and pitching moment (m P j ) totals are 
then computed by summing the local characteristics. j 



( 8 ) 

(9) 

( 10 ) 
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Figure 
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Mach Cone Correction Geometry 
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The strip rolling moment about the root chord is calculated from the airfoil 
lift and strip spanwise position, y : 


The total characteristics over the entire fin are then computed by summing 
the airfoil-strip characteristics in a spanwise direction. 


n 

■ L \ 

i = i 
j=i 

\ 

j=i 

■ £ \ 


The forces and moments computed in equations 4 through 15 are actually 
divided by the ambient dynamic pressure (q^) and strip width (Ay) . Thus, the 
associated coefficients are: 


2Fj • Ay 


(2 fins) 


2F d * Ay 


(2 fins) 


2M p ’ Ay 


(2 fins) 


4M r * Ay 


(4 fins) 


The center-of-pressure coordinates are calculated from the respective 
moments and the appropriate lift forces. 
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1 


X = 


J?S 

Fj cos a 


( 20 ) 


M 


Y = 


( 21 ) 


To obtain the linear slope coefficients near a = 0 at a given Mach number, 
the values of C L , C m , and Cj are calculated at a = 1° and are taken as the slope 
values CL a , C ma , and Ci g (all per degree). The slope values per radian are 
obtained by the proper conversion of radians to degrees. Thus: 


rad. 

C 

m q 

rad. 

C i 

h 

rad. 


180 

77 

( C L 

L . ,) 

(22) 

180 

(c 

i \ 

(23) 

77 

V m 

la = l) 


180 

77 

( c ll 


(24) 


The zero angle-of-attack values of C DW , X , and Y are simply their respective 
values at a = 0. 

CONFIGURATION AND GEOMETRY CONSIDERATIONS 

The user indicates the shape of the fin configuration by specifying the root 
chord length (C r ) , root airfoil leading and trailing edge region chord lengths 
(l Lr and l Tr ) , leading and trailing edge sweep angles (r L and r T ), region 
boundary sweep angles (r { and r 2 ), and the fin semispan length (S) (Figure 3). 
The program then computes the distance from the reference axis to the airfoil- 
strip leading edge in the streamwise direction (X L ) , airfoil chord length (C), and 
the airfoil region chord lengths ( 1 L andl T ) as functions of the spanwise coordi- 
nate of the airfoil strip. 


y tanP L cos 


a 


(25) 


C - y (tan P T - tan T L ) + C f 


(26) 
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SECTION A-A TYPICAL AIRFOIL 
Figure 3. Fin and Airfoil Geometry 
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(tanTj - tanT L ) + 1 L 

f 

(27) 

(tanr x - tariTj ) + 1 T 

(28) 


The width of the airfoil strip is determined from the number of strips (n) 
desired by the user. 


Ay - — (29) 

n 

The airfoil strip cros section is assumed to be symmetrical about the chord 
line. The most general crossection configuration is the modified double wedge 
(Figure la). The other three types (Figures lb, lc, and Id) are special cases of 
this basic shape. Figure 4 shows in detail the general crossection, which has 
six separate flow regions. As can be seen in this figure, the surface region 
lengths are determined by specifying the leading- and trailing-edge wedge half 


angles in the streamwise direction, 

£ L and £ T . 


1 1 ~ *4 

1 L 

(30) 

cos £ L 


1 3 = 1 . 


(31) 

J 6 

cos £ T 


* 2 - *5 

H 

•■H 

1 

^■4 

1 

O 

It 

(32) 


The program computes the local-surface inclination angles ( 75 J as functions 
of a : 



- * 

(33a) 

^2 = 

~ a 

(33b) 

^3 = 

- £ T - a 

(33c) 

^4 = 

+ a 

(33d) 
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7) s - a 

(33e) 

^6 = " : T + a 

(33f) 

From these, the streamwise and normal components of the region lengths 
are calculated. 

d i - 1. cos 7J. (i - 1 - 6) 

(34) 

n i = li sin Tj i (i = 1 - 6) 


The distance of the forward local-region boundaries from the reference axis 
are then determined as : 

X P1 X P4 \ 

(35a) 

X P 2 = X L + d l 

(35b) 

X P3 “ X L ^ d l + d 2 

(35c) 

X P5 " X L + d 4 

(35d) 

X P6 _ X L + d 4 + d 5 

(35e) 


FIN-TIP MACH CONE CORRECTION 

The Mach-cone angle (/x) is a direct function of the Mach number. 


Ar ctan ( 4 ) = Arcsin (-jj-) 


(36) 


For an airfoil strip at a spanwise distance y , the intersection of the strip and 
the Mach cone is a distance l w from the leading edge. 


1» " (S - y) tanP L + tan(90 - /x) 


(37) 
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But, by trigonometric identities and the use of equation 37, this may be reduced 
to: 


l w “ (S - y) (tanf L + /?) 


(38) 


If 1„ is greater than or equal to the chord length, there is, of course, no actual 
intersection and no correction is necessary. If 1 w is such that 


C > l w > C - 1 T 


then the ratio of the portion of 1 T falling behind the Mach cone to 1 T is : 

C - 1„ 


r 3 " r 6 


(39) 


Also, 


r l = f 2 - r 4 = r 5 = 0 


(40) 


If 1 w is such that 


C - 1 T > l w > i L 


then the corresponding ratio for the middle regions is 


C - 1 T l w 


r 2 r 5 


(41) 


Also, 


r 3 r 6 - 1 


(42a) 


and, 


f l = r 4 = 0 


(42b) 


If 1 w is such that 


c - 1 T - 1„ >1.10 


then the ratio for the leading edge regions is 


1l - It 


r l = r 4 


(43) 
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Also, 

r 2 “ r 3 - r S - r 6 - 1.0 

These ratios are the values used in equations 4, 5, and 6 for calculating the local 
lift, drag, and center of pressure. In these equations, the pressure coefficient in 
the portion behind the Mach cone is taken as half the value of the C p calculated 
by equation 1. 

PROGRAM USAGE 

The input data required by FIN are: 

M o Initial Mach number 

'>M Mach number increment 

M max Maximum Mach number 

Maximum angle of attack 

Reference area 
Reference length 

C r Root chord length 

S Fin semispan 

1 , Root airfoil leading-edge region chord length 

" r 

1 T Root airfoil trailing-edge region chord length 

Leading -edge sweep angle 

Region boundary sweep angles 
- sweep angle 

... ’ n the streamwi.se direction 



max 

A 

L 


» A 


'!• rect ion 




k = 0 


Number of strips desired 

Compute with fin tip Mach cone correction 


' 1 Compute without Mach cone correction 

Lengths and areas must be input in any consistant units. Angles are all input in 
degrees. 

Four input cards are required for each run. The input data are arranged in 
the order shown on the indicated card and the cards must be in the indicated 
order: 

• Card 1: Any 72 Hollerith characters for run identification 

• Card 2: «!,, AM, M„„, a.,,, A, L, C f , S 

$ Card 3: 1^ *^t 1 ^ ' 2 * ^t* 

r r 

• Card 4: n > k 

Any number of four-card sets may be input. 

Cards 2 and 3 contain eight fields each having ten spaces. Variables may be 
placed anywhere in their respective field, but each value on cards 2 and 3 con- 
tain a decimal point. Card 4 contains two fields each having three spaces. The 
values on card 4 do not contain a decimal point; however, the number must al- 
ways end in the third space in the field. 

Appendix B shows a typical set of inputs. 

The user may exercise a number of calculation options by proper definition 
of certain input parameters. If M 0 = 0 , then the program sets: 


where 


T" tan r L 


However, if < 0, the program sets M 0 = l. This option provides the minimum 
allowable Mach number for the given configuration. If a m>x = 1, the program 
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outputs only the values of C La , C DWq , C^, C 1? , X, and Y as functions of Mach 
number. As indicated in the input list, the fin-tip Mach cone may or may not be 
included by setting k equal to zero or one respectively. 

For a modified double wedge fin, the input values are simply those for the 
particular configuration. For the special cases of the modified double wedge, 
however, the configuration data must be input according to the following relations: 

Double Wedge 


r, = input value 

= input value 

r _ p 

2 1 

r T = input value 

Modified Single Wedge 
F L = input value 
T j = input value 

r = r 

2 T 

F T = input value 


Single Wedge 
P L = input value 

r : = r x 

r 2 = r T 

■ T = input value 


C r = input value 
1 L = input value 

r 



S L = input value 
S T = input value 


C r = input value 
1 L = input value 

r 

1 T =° 
r 

S L = input value 

S T =0 

C r = input value 

!«. = C r 

r 

h =° 

r 

S L = input value 

’T =0 
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RESULTS AND CONCLUSIONS 


To date, FIN has been applied to the fins of the Aerobee 150 A, Aerobee 350, 
and Tomahawk sounding rockets. Of these, only the first tv/o have enough relia- 
able windtunnel data to allow meaningful comparison (references 5 and 6). Un- 
fortunately, the data for these vehicles did not include enough parameters to 
allow a comparison of all the calculated characteristics. However, since all the 
characteristics stem from the same basic C p distribution, it is reasonable to 
infer the degree of accuracy, to a first approximation, of all the calculated 
quantities from the comparison of a few key characteristics. The characteristics 
that were compared to windtunnel data are C L , C L , and X o . The comparison 
calculations are not included in FIN. 


Both the Aerobee 150A and Aerobee 350 data contained total vehicle and 
body-alone data for C L , C L ^ , and X 0 . By combining the windtunnel body-alone 
data with FIN tail data, total vehicle values may be calculated for comparison 
with the total vehicle windtunnel data. The total vehicle parameters are obtained 
from the contributions of separate portions of the vehicle. 


C L I c l) b + ( C l) B (T) + ( C l) T (B) 

C l = (c L ) ♦ (c L ) + (c, ) 

u \ L a/B N L a/B(T) \ L a/T(B) 


The contributions of the body in the presence of the tail, (C L ) B(T) , (C L ) 
and(X o ) B(T) , are determined by the method of reference 2, using charts ( 
4b and 14b. The contributions of the tail in the presence of the body, (C L ) T(B) , 
( C L a )t ( b ) ’ 311(1 ( X o ) t ( b ) » are calculated from the FIN values (C L ) T , (C La ) T , 
and (X o ) T by applying the tail-body interference factor technique from refer- 
ence 2. 


(^l)t( B) 

= k t,b, ( c l) t 

( C Ot ( B) 

= K t<b,(c l J 

w*., 

= (*.) T 


( 47 ) 
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where K T(B) is the interference factor given by equation 14 and chart 1 of 
reference 2. Since this factor is valid only near a = 0, equation 47 is valid only 
for small angles of attack. Figure 5 shows C L vs. a for the Aerobee 150A at 
Mach 3.01. Figures 6 and 7 show C L vs. M for the Aerobee 150A and Aerobee 
350 respectively. Figures 8 and 9 show the respective curves of X o vs. M. 

As the figures show, the values calculated using FIN output are well within ^ 

10% of all the windtunnel values at low angles of attack between Mach 2 and 7. 

This represents adequate prediction for use in first-look trajectory analyses or 
for determining requirements for windtunnel tests on a new or proposed sounding 
rocket vehicle. 
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Figure 6. Lift Coefficient Slope as a Function of Mach Number, Aerobee 150A 
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Figure 9. Center of Pressure from Nose as a Function of Mach Number, Aerobee 350 
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dec -j a oui^oj nh 000000000000000 o o o o <* 
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$JOB 1091P003 405 J SNOW BARROWMAN 
$EXECUTE I B JOB 

IBJOB GO, SOURCE, NOMAP 

IBFTC FIN M94, XR7 » DECK 

FIN 

CALCULATION OF, 

CL X ALPHA 
CM X ALPHA 

CL X DELTA (ROLLING MOMENT COEFF.) 

COW X ALPHA 

CENTER OF PRESSURE X ALPHA 
LIFT CURVE SLOPE 
PITCHING MOMENT CURVE SLOPE 
ROLLING MOMENT CURVE SLOPE 
WAVE. ORAG CURVE SLOPE 

OF A THREE DIMENSIONAL FIN. PROGRAM USES BUSEMANNS SECOND OROER 
THEORY WITH THE THIRD ORDER TERMS RETAINED. TWO DIMENSIONAL VALUES 
ARE CALCULATED ALONG STREAMWISE STRIPS AND ARE THEN SUMMED IN A 
SPANWISE DIRECTION. THE EFFECT OF THE FIN TIP MACH CONE IS 
CONSIDERED BY DECREASING THE LOCAL PRESSURE COEFFICIENT BY A 
FACTOR OF 1/2 WHEN THE POINT IS WITHIN THE MACH CONE. 

GOOO FOR WEDGE, DOUBLE WEDGE, OR MODIFIED DOUBLE WEDGE AIRFOILS 
WITH SYMMETRY ABOUT THE CHORD LINE. 

DIMENSION ETA(6),ETAR(6),CD0( 100 ) ,CLA ( 100 ) ,TLAM ( 4 ) ,L AM ( 4 ) ,N ( 6 ) 
DIMENSION XF (6) ,CLD( 100) ,CMA< 100 ) , CPXX ( 100 ) ,CPSS( 100) 

DIMENSION PL(6),D(6),R(6),FL(6),FD(6),TITLE(12),M(100) 

REAL M,N,LAM,LW,M0,K1,K2,K3,LL,LT,LIFT »MMAX »MSQ , LAR ,MU 

REAL LLR , LTR , MO »MMOMT , LMOM 

DTR=3. 1415927/180.0 

READ (5,2) (TITLEU ), 1-1,12) 

FORMAT (12A6) 


C 


READ (5,3) MO,DELM,MMAX, AMAX, AREA, REFL »CR, SPAN 
READ (5,3) LLR, LTR, (LAM( I ) , 1=1,4) ,ZETAL*ZETAT 
FORMAT (BE 10.4) 

READ (5,4) NS, 12 
FORMAT (213) 

WRITE (6,5) 

FORMAT (1H1) 

WRITE (6,2) (TITLEU ) ,1 = 1,12) 

F fN ANGLE VALUES 
DO 6 1=1,4 
LAR=LAM(I )*DTR 
TLAM ( I )=SIN(LAR)/COS(LAR) 

ZETALR=ZETAL*DTR 
ZETATR=ZETAT*DTR 
CZL=COS(ZETALR ) 

CZT=COS(ZETATR) 

WRITE (6,7) ZETAL ,ZETAT 

FORMAT ( 1 1H0 ZETA L * F6.3»5H DEG. , 5X, 9HZETA T = F6.3,5H DEG.) 
WRITE (6,e) AREA, REFL, NS 

FORMAT ( 1 3H0 AREF = F8.2,9H LREF = F8.3,UH STRIPS = 13) 

STRIP WIDTH 


XNS=NS 
DS=SPAN/XNS 
XNS=2.*DS/AREA 
INITIAL MACH NUMBER 
J=1 

IF (MO.NE.O.) G0T010 
CMU*CR/SPAN-TLAM( 1) 

IF (CMU.LE.O. ) G0T09 
M( J)=S0KT(CMU*CMU+1. ) 
GO TO 12 
M ( J ) = 1 .0 
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GO TO 12 
M( J)=MO 
GO TO 12 
MACH NUMBEK 
JsJ+1 

M< J)=M( J-D + DELM 

IF (M( J) .GT.MMAX ) G0T043 

MSO=M( J)*M( J) 

BETA=SOKT(MSU-l. ) 

DEM=-2.*MSg+4./3. 

DEN=l./bETA**7 

K1=2./BETA 

K2=( 1 . 2*MSg**2-2.*BETA**2 )/BETA**4 

K3=( .4*MSg**4-( 10.88/6. ) *MSU** 3+4. *MSQ*MSQ+DEM ) ♦DEN 
83= I .04*MSg**4-. 32#MS0**3+. 4*MSg*MSQ) *DEN 
MACH ANGLE 

CMU=(TLAM( 1 J+BETA) 

MU=SPAN*CMU 
IF (AMAX.EQ.l.) GOTO 16 
WRITE (6,13) 

FORMAT (1H0) 

WRITE (6,14) M ( J ) 

FORMAT ( 5H0M = Fb.3) 

WRITE (6,15) 

FORMAT (58H ALPHA CL CM CLM COW 

1CPS) 

;l VS. alpha CURVE 
N AM AX = AM AX + 1 • 

GO 42 K= 1 ,NAMAX 
ALPHA=K-1 
ALPHAR=ALPHA#l)TR 
TL I FT =0 . 

TDRAG=0. 

TMM0M=0 • 

TLMUM=0 • 

S=-OS 

S=$«-DS 

C=S*(TLAM(4)-TLAM( 1) )+CR 

LL=S#(TLAM(2)-TLAM( 1) )+LLR 

LT=$*(TLAM(4)-TLAM{ 3) )+LTR 

Sl=LL/CZL 

ST-LT/CZT 

SM=C-LL-LT 

XL = S*TLAM'( 1 )*COS( ALPHAR ) 

ETA< 1)=ZETAL-ALPHA 
ETA( 2 ) “-ALPHA 
ETA(3)=-ZETAT -ALPHA 
ETA(4)*ZETAL+ALPHA 
ETA ( 5 ) =ALPHA 
ETA(6)=-ZETAT+ALPHA 
DU 18 1=1,6 
ETAR ( I ) =ETA ( I ) *DTR 
D ( 1 ) = SL*COS ( E TAR ( 1 ) ) 

D(2)=SM*C0S(ETAR (2) ) 

Dl3)=ST*C0S(ETAR(3) ) 

D ( 4 ) = SL*COS ( ET AR (4 ) ) 

D ( 5 ) = SM*COS ( E TAR ( 5 ) ) 

D 1 6 ) = $ T*COS ( ET AR ( 6 ) ) 

N( 1) = SL*SIN(ETAR( 1) ) 

N ( 2 ) = SM*SIN ( ETAR ( 2 ) ) 

M3 ) “ST* SIN I ETAR ( 3) ) 

N(4)=SL*SIN(ETAR(4) ) 

M 5 ) “SM*SIN (ETAR ( 5) ) 

N(6)=ST*SIN(ETAR(6) ) 

DO 21 1=1,6 

SAVE=K1*ETAR(| ) +R2*ETAR ( I )**2*K3*ETAR( !)»*3 


IF II.GE.A) G0T020 
IF lETAR(l).GT.O.O) G0T019 
R<I)=SAVE 
GO TO 21 

19 R(I)»SAVE-B3*ETAR( 1)**3 
GO TO 21 

20 R ( I )*SAVE-B3*ETAR ( A) **3 

21 CONTINUE 

IF (I2.NE.0) G0T028 
C FIN TIP MACH CONE CORRECTION 
LW*MU-S*CMU 
IF (LW.GE.C) G0T026 
SL*C-LT 

IF (LW.GE.SL) G0TQ2A 
IF (LM.GE.LL) G0T023 
PL<1)=(LL-LW)/LL 
PL (A )*PL II) 

DO 22 I =2*3 
PL ( 1 + 3) * 1 • 

22 PLI I )*1» 

GO TO 30 

23 PL(2)*<C-LT-LW)/SM 
PL(5)*PL<2> 

PL(3)*1. 

PL (6 )*PL I 3) 

PL( 1)=0. 

PL ( A )*0. 

GO TO 30 

2A PL ( 3 )* (C-LW ) / LT 

PL ( 6 ) *PL ( 3 ) 

DO 25 1 = 1,2 
PL ( I )=0. 

25 PL ( 1 +3 ) =0 • 

GO TO 30 

26 DO 27 I * 1 »6 

27 PHI ) *0. 

GO TO 30 

28 00 29 1 * 1,6 

29 PL (I) *0.0 

C STRIP COEFFICIENTS 

30 00 31 I « 1 ♦ 3 

31 FLU )*-K ( I l*DII)*l l.-.5*PL(I )) 

00 32 I *A ,6 

32 FLIl)>Rin*0U)*ll.-.5*Pim) 

00 33 1 = 1*6 

33 FDII)*RII)*Nin*ll.-.5«PLIl)) 

00 3A 1*1, A, 3 

3A XP(I)*XL 

00 35 1*2, 5, 3 

35 XP(I)*XL+D(1-1> 

DO 36 1*3,6, 3 

36 XPm*XL+0(I-l)+0l 1-2) 

MMOMT *0. 

00 37 1*1,6 

CW*.5*D( I )*< l.-PLI I ) ♦• 5*PL < I ) **2*XP ( I ) *1 2. -PL II ) )/0(I ) ) /I l.-.5*PL< 
11 )) 

37 MMOMT *MMOMT-CW*FL ( I ) 

LIFT *0*0 

DRAG*0.0 
00 38 1*1,6 
LIFT*L1FT+FL( I ) 

38 DRAG-ORAG+FOtl ) 

LMOM«S*LIFT 

TLIFT-TLIFTaLIFT 

TORAG-TDKAG+DRAG 

TMMOM-TMMOM+MMOMT 


ill I 


P 

i 


TLHOM=TLMOH+LMOM 
IF (S.LT.SPAN) G0T017 
C TOTAL COLFF I C I ENTS 
CL=TL IFT*XNS 
CD=TDRAG*XNS 
CM=TMMOM*XNS/REFL 
CLM=2.*TLMOM*XNS/REFL 
CPX=TMMOM/(TLIFT*COS( ALPHAR ) ) 
CPS=TLMUM/TL1FT 
IF (K.EU.l) CDO(J»=CO 


IF (K.EU.2) G0T039 
GO TO 40 

C LIFT CURVE SLOPE 

39 CLAI J )=CL/DTR 
CLl)(J)=CLM/DTR 
CMA< J) = CM/OTR 
CPXXI J)=CPX 

r dqqi it = rp c 

40 IF (AMAX.EQ.l.) G0T042 

WRITE (6,41) ALPHA,CL,CM,CLM,CD,CPX,CPS 

41 FORMAT ( 1H F9 . 3, 4F8 . 4, 2F9 . 3 ) 

42 CONTINUE 
GO TO 11 

43 J=J-1 

IF (AMAX.EO. 1. ) G0T045 


WRITE (6,44) 







FORMAT (63H1 
CPS) 

GO TO 47 
WRITE (6,46) 

M 

CLA/RAD. 

CMA/RAD. 

CLD/RAD. 

COW 

CPX 

FORMAT ( 6 3H0 
CPS) 

M 

CLA/RAD. 

CMA/RAD* 

CLO/RAO. 

COW 

CPX 


47 WRITE (6,48) ( M ( I ) , CL A ( I ) , CMA ( I ) , CLDII) ,CDO ( I ) ,CPXX (I ) ,CPSS ( I) , I * 1 
1, J) 

48 FORMAT ( 1H F7 . 3, F9 . 3, 2F 10 . 4 , 3F9. 4 ) 

GO TO 1 

END 


* 


( 

i 

1 

( 


30 


APPENDIX B 
SAMPLE CASE 

Figure B-l shows the Tomahawk sounding rocket fin. For this run, the in- 
put data were chosen as follows. 



= 3.0 

AM 

= 1.0 


= 7.0 


= 5.0 degrees 

A 

= 63.6 in. 2 

L 

= 9.0 in. 

C , 

= 21.9 in. 

S 

= 13.8 in. 

l Lr 

= 3.76 in. 

*Tr 

= 0. in. 

r L 

= 45 degrees 

r i 

= 38.23 degrees 

r a 

= 0. degrees 

r T 

= 0. degrees 

Cl 

= 2.475 degrees 


= 0. degrees 

n 

= 300 

k 

= 1 


Figure B-2 shows a listing of the input cards. 
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The output resulting from these inputs is shown below. 

TOMAHAWK FIN INPUTS 

2ETA L = 2.475 DEG • ZETA T = 0 . DEG. 

AREF = 63.60 LREF = 9.000 STRIPS = 300 f 


H = 3.000 


ALPHA 

CL CM 

CLM 

CDW 

CPX 

CPS 

0 . 

0.0000 -0.0000 

0.0000 

0.0029 

-7.695 

5.834 

1.000 

0.1655 -0.2519 

0.2147 

0.0058 

-13.702 

5.839 

2.000 

0.3313 -0.5040 

0.4290 

0.0146 

-13.702 

5.839 

3.000 

0.4976 -0.7566 

0.6457 

0.0293 

-13.702 

5.839 

4.000 

0.6648 -1.0098 

0.8627 

0.0499 

-13.703 

5.839 

5.000 

0.8332 -1.2638 

1.0812 

0.0767 

-13.703 

5.839 


t-i = 4.000 


ALPHA 

CL CM 

CLM 

CDW 

CPX 

CPS 

0 . 

O.UOOO -0.0000 

0.0000 

0.0021 

-5.407 

5.340 

1.000 

0.1219 -0.1847 

0. 1580 

0.0043 

-13.641 

5.836 

2.000 

0.2442 -0.3699 

0.3167 

0.0108 

-13.642 

5.836 

3.000 

0.3674 -0.5562 

0.4765 

0.0217 

-13.642 

5.836 

4.000 

0.4921 -0.7442 

0.6382 

0.0372 

-13.644 

5.836 

5.000 

0.6187 -0.9345 

0.8024 

0.0572 

-13.645 

5.836 


r*. = 5.000 


ALPHA 

CL 

Ci*i 

CLM 

CDW 

CPX 

CPS 

0 . 

0.0000 

-0.0000 

0.0000 

0.0017 

-6.606 

5.268 

1.000 

0.0972 

-0.1466 

0.1259 

0.0035 

-13.579 

5.832 

2.000 

0.1950 

-0.2940 

0.2527 

0.0087 

-13.580 

5.832 

3.000 

0.2941 

-0.4431 

0.3811 

0.0175 

-13.581 

5.832 

4.000 

0*3951 

-0.5949 

C.5121 

0.0300 

-13.583 

5.832 

5.000 

0.4987 

-0.7500 

0.6464 

0.0463 

-) 3.586 

5.833 


n = 6.000 


ALPHA 

CL CM 

C L H 

CDW 

CPX 

CPS 

0. 

0.0000 -0.0000 

0.0000 

0.0015 

-7.423 

5.767 

1.000 

0.0812 -0.1219 

0.1051 

0.0029 

-13.515 

5.828 

2.000 

0.1631 -0.2449 

0.2113 

0.0073 

-13.516 

5.828 

3.000 

0 • 24o 7 -0.3701 

0.3196 

0.0148 

-13.519 

5.829 

4.000 

0.3328 -0.4988 

0.4M0 

0.0254 

-13.523 

5.829 

5.000 

0.4221 -0.6320 

0.5468 

u.C3<H 

-13.529 

5.829 


i = 7.000 


ALPHA 

CL CM 

CLf*i 

COW 

CPX 

CPS 

to 

0. 

0.0000 -0.0000 

0.0000 

0.0013 

-7.159 

5.629 

* 

1.000 

0.0699 -0.1045 

0.0905 

0.0026 

-13.449 

5.824 


2.000 

0.1408 -0.2104 

0. 1823 

0.0064 

-13.452 

5.825 


3 .000 

u • 2 1 37 -0.3191 

0.2766 

0.0129 

-13.457 

5.825 


4 . 00 0 

0.2894 -0.4319 

0.3747 

U .0222 

-13.^64 

5.825 


5.000 

0.^691 -0.5505 

0.4779 

0.0346 

-13.472 

5.826 



M 

cla/kad. 

CMa/KAD. 

CLD/RAD. 

Ci>H CPX 

CP i 

5.000 

9.482 

-14.4327 

12. 3032 

0.0029 -13.7018 

5.8392 

4. GUO 

6.982 

-10.5806 

9.0540 

0.0021 -13.6413 

5 .8356 

5.000 

5.567 

-8.3980 

7.2149 

0.0017 - 1 3 .5786 

5.8 52 * 

6*000 

4.651 

-6.9823 

6.0232 

0.0015 -13.5145 

5.6282 

7.000 

4.U07 

-5.9666 

5. 1859 

0.0013 -13.4494 

5.8/44 
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